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Abstract 

> 

Since the complexity of the practical environment, many distributed networked 
0^ , systems can not be illustrated with the integer-order dynamics and only be de- 

scribed as the fractional-order dynamics. Suppose multi-agent systems will show 
the individual diversity with difference agents, where the heterogeneous (integer- 
• order and fractional-order) dynamics are used to illustrate the agent systems and 

compose integer-fractional compounded-order systems. Applying Laplace transform 
and frequency domain theory of the fractional-order operator, consensus of delayed 
multi-agent systems with directed weighted topologies is studied. Since integer-order 
model is the special case of fractional-order model, the results in this paper can be 
extend to the systems with integer-order models. Finally, numerical examples are 
used to verify our results. 
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1 Introduction 



Nowadays, the past three decades have witnessed significant progress on frac- 
tional calculus, because the applications of fractional calculus are found in 
more and more scientific fields, covering mechanics, physics, engineering, infor- 
matics, and materials. The list of such applications is long, for instance, it in- 
cludes viscoelasticity, colored noise, dielectric polarization, electrode-electrolyte 
polarization, electromagnetic waves, control engineering and so onjTE '■'■> 1 (i 71IH1. 
In fact, real word processes generally or most likely are fractional-order systems [TP] . 
Furthermore, fractional order controllers have so far been implemented to en- 
hance the robustness and the performance of the closed loop control systems, 
the stability problem of fractional-order systems has been investigated both 
from an algebraic and an analytic point of view[ 



On the other hand, as one of the most basic problem of coordination control 
for networked systems, consensus of multi-agent systems has been widespread 
concerned as an important research topic in the field of systems control. Con- 
sensus of multi-agent systems means that several distributed agents achieve 
the same state or output through local mutual coupling effect among the in- 
dividuals, where the centralized control is not used. Based on the computer 
model proposed by Reynolds which imitates animals' flocking [9], Vicsek et 
al. [10] firstly proposed a non-equilibrium multi-agent system model from the 
point of view of statistical mechanics, simulation shows that all individuals 
in the system can run in accordance with the same direction under certain 
conditions. Recently, since the wide application of multi-agent systems in var- 
ious fields, many scholars have devoted themselves to study the consensus 
of multi-agent systems [Ilfl2|ll3|ll4|ll5fl6fl7j . When the information transfer- 
ring, communication delays will occur in networked control system and affect 
the features of the system. The effects of communication delays on move- 
ment consensus of multi-agent systems have been concerned by many scholars 
|14|ll5|ll6|fTTj . and the stability of delay system also become a hot topic in the 
multiple agents field. 

The important results of the above literatures pay attention to the consensus 
problem of integer-order multi-agent systems. In the complex environment, 
many dynamic characteristics of natural phenomena can not be described in 
the form integer-order equation, but only be described in the dynamics of 
fractional-order (non-integer order) behavior, for example: flocking movement 
and food searching by means of the individual secretions and microbial, sub- 
marine underwater robots in the bottom of the sea with a large number of 
microorganisms and viscous substances, unmanned aerial vehicles running in 
the complex space environment [IS] . Cao and Ren [191120] studied distribu- 
tion coordination of fractional-order multi-agent systems firstly, and gave the 
relationship between the number of individuals and the fractional order in 
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the stable multi-agent systems. However, to the best of authors' knowledge, 
there are few researches done on the coordination control of fractional-order 
multi-agent systems with communication delays. 

In this paper, we suppose agents work in the complex environment, the het- 
erogeneous dynamics with fractional-order and integer-order is presented for 
multi-agent systems. The main innovation of this paper lies in the study on 
consensus of compound-order (fractional-order and integer-order) distributed 
multi-agent systems with different time delays. This paper is organized as 
follows. In Section 2, some necessary definitions and notations are given on 
fractional calculus. In Section 3, the compound-order dynamics of the frac- 
tional order system and the integer order system is presented. The consensus 
of integer-fractional compounded-order multi-agent systems with communi- 
cation delays is studied in Section 4. The corresponding simulation results 
are provided in Section 5 to demonstrate the effectiveness of the proposed 
conditions. Finally, the conclusions are drawn in Section 6. 



2 Fractional Calculus 

Fractional calculus plays an important role in modern science. There are 
mainly two widely used fractional operators: Caputo and Riemann-Liouville 
(R-L) fractional operatorspQ. In physical systems, Caputo fractional operator 
is more practical than R-L fractional operator because R-L fractional oper- 
ator has initial value problems. Therefore, in this paper we will use Caputo 
fractional operator to model the system dynamics and analyze the stability of 
the proposed coordination algorithms. Generally, Caputo fractional operator 
includes Caputo integral and Caputo derivative. Caputo integral is defined as 

where the integral order p £ (0, 1], T(.) is the Gamma function, and a is an 
arbitrary real number. Based on the Caputo integral, for a nonnegative real 
number a, Caputo derivative is defined as 

C a Dff(t) = c a D?[^f{t)], (1) 

where p = [a] + 1 — a G (0, 1] and [a] is the integral part of a. If a is an 
integer, then p = 1 and the Caputo derivative is equivalent to the integer- 
order derivative. In this paper, a simple notation f^(t) is used to replace 
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^Dff(t). Let £() denote the Laplace transform of a function, the Laplace 
transform of Caputo derivative is shown as 

£(/<«>(*)) = s<*F(s) - S'a-V^HO). (2) 

fc=i 

where F(s) = £(/(£)) = f^e~~ st f(t)dt is the Laplace transform of function 
f(t), /«(0) = lim^o- / (fe) (0 and /(°>(0) = /(0) = lim^ - 



3 Problem statement 



Assume that multi-agent systems consist of n autonomous agents, connected 
relations among the agents constitute a network topology Q. Let Q = {V, E, A} 
represent a directed weighted graph, in which V = {vi, v 2 , v n } represents 
a collection of n nodes, and its set of edges is E C V x V. The node indexes 
belong to a finite index set I = {1, 2, n}, adjacency matrix A = [ay] G _R nxn 
with weighted adjacency elements a»j > 0. An edge of the weighted diagraph 
Q is denoted by = G £7. We assume that the adjacency element 

dij > when G otherwise, a^- = 0. The set of neighbors of node i is 
denoted by Ni = {j e I : > 0}. 

Let Q be a weighted digraph without self- loops, i.e., = 0, and matrix 
D = diag{di, c?2j •••> ^n} be the diagonal matrix with the diagonal elements 
di = Dj=i a ij representing the sum of the elements in the i-th row of matrix 
A. The Laplacian matrix of the weighted digraph Q is defined as L = D — A. 
For two nodes % and k, there is subscript set {k\, k 2 , ...ki) satisfying > 0, 
o-fciAi2 > 0, a^k > 0, then there is a directed linked path between node i and 
node k which is used for the information transmission, also we can say node 
i can receive the information from node k. If node i can find a path to reach 
any node of the graph, then node i is globally reachable from every other node 
in the digraph. For any two nodes in the graph, there are at least one directed 
linked path, then G is strongly connected. 

Lemma 1[I2]- is a simple eigenvalue of Laplacian matrix L, and X = 
C[l, 1, 1] T is corresponding right eigenvector, i.e., LX = 0, if and only if 
the digraph Q = (V, E, A) has a globally reachable node. 

Assume that there are individual differences in the complex environment of 
multi-agent systems; there are two groups of multi-agent systems with integer- 
order dynamics and fractional-order dynamics. The compounded-order dy- 
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namical equations are described as follows: 



Xi(t) = Ui(t),i = 1, ...,m, 
x\ a \t) = ui(t), I — m + 1, ...,n, 



(3) 



where Xi(t) G R and Ui(t) G R represent the i-th agent's state and control 
input respectively, Xi(t) represents the first-order derivative for the state Xi(t), 
x\ a ^ represents the a order Caputo derivative, and a G (0,1]. Assume the 
following control protocols are used in multi-agent systems: 

Ui(t) = -7 a ik[%i(t) - x k (t)],i = 1, ...,m,m+ 1, ...,n. (4) 



where represents the (i, k) elements of adjacency matrix A, 7 > is control 
gain, iVj represents the neighbors collection of the i-th agent. 

This article assumes that there are communication delays in the dynamical sys- 
tems, and consensus of the integer-fractional-compounded-order agent systems 
with communication delays is studied. Under the influence of communication 
delays, we can get the following algorithm: 

Xi(t) =Ui(t-Ti),i = l,...,m, 
x\ a \t) = u t (t -Ti),l = m + 1, ...,n, 

where > is the communication delay of agent i. Through a simple change 
we can get 

[*!(*), xjrY = -i^? =l L t x T {t - n ), ^U+.UX 1 \t - n )] T (6) 



where Lj = E(i)L, E(i) represents matrix whose element of is 1 and 
the rest are 0s, X^t) = [x^t), x m (t)], X 2 (t) = [x m+1 (t), x n (t)], X(t) = 
[xi(t), x 2 (t), x n (t)], and L = S" =1 Lj. Suppose that for any initial value of 
the system, the states of autonomous agents meet lim t ^. 00 (x i (t) — Xk(t)) = 0, 
for i,k G /, then we call multi-agent systems asymptotically reach consensus. 
In this paper, we apply Laplace transform and frequency domain theory of the 
fractional-order operator to study the consensus of delayed compounded-order 
multi-agent systems with directed weighted topologies. 
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4 Coordination control for compounded-order multi-agent systems 
with communication-delays 



Theorem 1 Suppose that multi-agent systems are composed of n indepen- 
dent agents whose connected network topology is directed and has a globally 
reachable node. Then, compounded-order multi-agent system ([6]) with time 
delays can asymptotically reach consensus, if 

Tl < 2(2 7 J)V" (7) 



where d = max{<ij, i = 1, n}, di = S^ =1 a ifc . 



Proof. Applying Laplace transformation to system 
equation can be gotten 



), the characteristic 



det( 



Sir, 



\ 



J-n—rn I 



+ 1 E{s)L) = ^ 



where I m represents a unit matrix with m-dimensions, E(s) = diag{e _TlS , e~ TnS }. 
is a single eigenvalue of the Laplacian matrix L because the system has a 
globally reachable node from Lemma 1. Due to a > 0, the characteristic equa- 
tion has a characteristic root s = 0. 



When s ^ 0, let 



F( S ) = (J n + 7 



V 



E(s)L), 



the characteristic equation (jSJ) is equivalent to F(s) = 0. Nextly, we will prove 
that all zero points of F(s) = have negative real parts. Let 



/ 



G(s) = 7 



V 



s- Q /„ 



E(s)L, 



according to the generalized Nyquist criterion |21j . if for s = ju>, where j is 
complex number unit, point —1 + jO is not surrounded by the Nyquist curve 
of G(juys eigenvalues, then all zero points of F(s) have negative real parts. 
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Let s = ju>, we can get 

-lp-JTT/2 J \ 

w -a e -W2 /n _ m ^ 



G(jw) = 7 



(9) 



In the following proof, Gerschgorin's disc theorem will be applied to estimate 
the eigenvalues X(G(ju)) of the matrix G(juo). According to the Gerschgorin's 
disc theorem, we have 

A(G(j'w))eU G « ( 10 ) 
iei 



where for i — 1, m, 



G t = {( e C, |C - w- 1 7 d i e-^ WTi+ ' r / 2 ) 
and for i — m + 1, n, 



= {( G G, |C - w- Q 7^e-^ WTi+a7r / 2 )| 

where = Sfe=i Qife- After simple sorting it can be 

Gi = {( e C, |C - w-Vie- J ' (wTi+,r/2) | < = 1, 

and 

Gi = {C e C, |C - w- Q 7^e-' J '^ +a ^| < co~ a -fdi},i = m + 1, ...,n, 



When the Nyquist curve of the origin of the disc G, changing, the disc changes 
along with it. Next, we will prove that point —a + jO(a > 1) is not in every 
disc Gi. 

For i — 1, ...,m, the changes of the disc Gj is to be studied in the following. 
The origin of the disc G« is co" 1 'jdie~^ UJTi+7T ^ 2 \ the radius of the disc Gj is 
w _1 7 cij. Let 
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A = \-a + jO- urV;e- j( " ri+7r/2) | 2 - (w"^) 2 , 
there is 

A = a(a + 2u~ lr ydi cos(wTj + 7r/2)). 
When w c Tj + 7r/2 = 7r, there is cos(w c Tj + 7r/2) = —1. We can get 

A > a(a - 2u~ 1 ^d i ) ) 
where 

2a; 6 7 V* = 2(7r/(2r,))" 1 7^. 
According to the conditions of the theorem 

2(7r/(2r i ))- 1 1 d i < 2((ir/(2r i ))- a 1 d< 1, 
from the hypothesis a > 1, we can get 

A > 0. 

Then, it gets, for i — 1, ...,m, 

| - a + j'O - urV*e- j( " ri+7r/2) | > w-^. 
By means of the same deduction, we can get, for i — m + 1, n, 

\-a + j0- uj- a ^die- j ^ Ti+a7T/2) \ > uj-^di. 

When a > 1, the point —a + jO is not in disc Gj. Thus the point —1 + jO is 
not surrounded by curves of eigenvalue X(G(juj)) of matrix G(juj). Therefore, 
all zero points of F(s) = have negative real parts. 

Due to the equilibrium point of the system meeting LX* = 0, then X* = 
C[l, 1] T (where C is a constant) is the eigenvector with the corresponding 
eigenvalue of the Laplacian matrix L. Therefore, lim^oo Xi(t) = C, and the 
system asymptotically reaches consensus. 

Corollary 1. Suppose multi- agent systems are composed of n independent 
agents, whose connection network topology is directed and symmetrical, and 
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there is a global reachable node. Then compounded-order multi-agent system 
OH]) with time delays can asymptotically reach consensus, if 



Ti < 



7T 



2( 7PL )iA 



(11) 



where pi represents the spectral radius of matrix L with pi = max{|Aj|,i = 
1, ...,n} and Aj is the eigenvalue of the Laplacian matrix L. 

Proof. According to theorem 1, the characteristic equation of the system is 



det( 



Sin 



\ 



+ 1 E{s)L) = 0. 



Since the Laplacian matrix L is symmetrical, there is orthogonal matrix P 
satisfying L = PAP -1 , where A = diag{\i, X n }. Because there is a global 
reachable point, we can know Rank(L) — n — 1 and is a single eigenvalue of 
matrix from Lemma 1. Therefore, the characteristic equation has a root s = 0. 

When s ^ 0, let P(s) and G(s) be same as the proof of Theorem 1, and let 



Ms) 



S~ l I 



s~ a I 



diag{e ns ,i = 1, ...,n}. 



(12) 



Let s = jco, we can get 
G(ju) = H(ju>)>yL, 



(13) 



where 



H(ju) = di&g{Hi(ju)),i = 1, n} 



-jair/2 j 



diag{e juJTi ,i = 1, ...,n}. 



Let M = diag{Mj, % = 1, ...n} where M» = it/(2t^ (for % = 1, m) and M l = 
((2-a)7r/(2T l )) a (for / = m + l,...,n). Matrix MH(juj) = diag{MiHi(ju), i = 
1, n} is a diagonal matrix where the Nyquist curve of its diagonal elements 
passes over point —1+jO. Suppose \(G(jcj)) is the eigenvalue of matrix G(jcu), 
we have 
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X(G(ju)) = \(MH(ju)'yM~ 1 / 2 LM 1 / 2 ) 

G p( 7 M- 1 / 2 LM 1 / 2 ) x Co(QU{MiHi(ju),i= l,...,n}), 

Where p() represents the spectral radius of matrix, Co(£) represents the convex 
hull of £. Because of MiHi(Ju) will passes over point —1 + jO, point —1 + jO 
is included in convex hull Co(0 U {MiHi(ju),i = 1, ...,n}). Since 

p( 7 M" 1 / 2 LM 1 / 2 ) = p( 7 M- 1 / 2 AM 1 / 2 ) 

= max{| 7 Mj _1 Aj|, z = 1, n}, 

according to the hypothesis condition 7* < 2( - 7p ^i /Q , we can get 

p( 7 M~ 1 / 2 LM 1 / 2 ) < 1. 

Therefore, point -1+jO is not included in p(jM- l / 2 LM 1 / 2 )xCo(0U{M i H i (juj), i = 
1, ...,n}). That is, point —1 + jO is not included in the Nyquist curve of the 
eigenvalue of G(ju). According to generalized Nyquist theorem [2TJ, the zero 
points of the characteristic equation have negative real parts. Therefore the 
multi- agent systems can asymptotically reach consensus, and lim^oo Xi(t) = 
C. 

Corollary 2. Suppose multi- agent systems are composed of n independent 
agents, whose connection network topology is directed and symmetrical, and 
there is a global reachable node. Then compounded-order multi-agent system 
([6]) with time delays can asymptotically reach consensus with a — 1, if 

7T, < 7v/(2\ max ), (14) 
where X max is the maximum eigenvalue of matrix L. 

Corollary 3. Suppose multi- agents are system composed of n independent 
agents, whose connection network topology is directed and symmetrical, and 
there is a global reachable node. Then compounded-order multi-agent system 
(jBJ) can asymptotically reach consensus when a = 1 and time delays T{ = r, if 

2 7 r < 7r/p, (15) 
where p represents Spectral radius of matrix L. 

Remark 1. The consensus result in Corollary 3 for 7 = 1 is in accord with 
that in pi]. 
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5 Examples Simulations 



Suppose the system is composed of four dynamical agents (Fig. 1) with two 
integer-order agent systems (agnetl and agent2) and two fractional-order 
agent systems (agent3 and agent4). The connection weights between indi- 
viduals are 021 = 0.7, = 0.8, = 0.9, au = 1. The order of the fractional 
multi-agent dynamics is a = 0.9, through the topology of the system,we can 
get the adjacency matrix 



1^ 
0.7 
0.9 

0.8 Oy 




Fig. 1. Network topology of the multi-agent systems. 



According to Theorem 1, we can get the relationship between the system con- 
trol gain and the upper bound of communication delays (Fig. 2). With the help 
of Fig. 2, we can select the control gain according to the communication delay 
of the system, or decide the upper bound of communication delays by means 
of the control gain of the system, to make the system meet the condition of 
reaching consensus. Suppose the communication delay is r = 0.6s, the system 
control gain should be selected as 7 < 1.19 from Fig. 2; suppose the system 
control gain 7 = 1, we can obtain that the upper bound of communication 
delays is 0.7s from Fig. 2. 
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12 



Fig. 2: Relationship between the control gain and the upper bound of 

communication delays. 




40 60 
time(s) 



100 



Fig. 3: Movement trajectories of the multi-agent systems with delay 0.6s. 

Assume the communication delay of multi-agent systems is 0.6s and the sys- 
tem control gain 7 = 1 in simulation, we set the expect objective at 0.5, 
consensus can be asymptotically reached (Fig. 3) through compounded-order 
coordination algorithm. 



Assume the communication delay of multi-agent systems is 0.7s and the system 
control gain 7 = 1 in simulation, consensus of compounded-order multi-agent 
systems can not be reached (Fig. 4). 
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0.8 




time(s) 

Fig. 4: Movement trajectories of the multi-agent systems with delay 0.7s. 
6 Conclusions 



This paper studies distributed coordination of integer-fractional-compounded- 
order multi-agent systems with communication delays. Consensus of multi- 
agent systems with directed network topology is studied through the stability 
theory of frequency domain, and the consensus conditions for compounded- 
order delayed multi-agent systems are presented. The relationship between 
the control gain of multi-agent systems and the upper bound of time delays 
is derived. Suppose the orders of the fractional dynamical systems are all 
1, the extended conclusion in this paper is the same with ordinary integer 
order system. In the following work, research of the robust stability of integer- 
fractional-compounded-order multi-agent systems will be carried out. 
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